Nature of the quantum phase transition to a spin-nematic phase 
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It is shown that the quantum phase transition in metallic non-s-wave ferromagnets, or spin ne- 
matics, is generically of first order. This is due to a coupling of the order parameter to soft electronic 
modes that play a role analogous to that of the electromagnetic vector potential in a superconductor, 
which leads to a fluctuation-induced first-order transition. A generalized mean-field theory for the 
p-wave case is constructed that explicitly shows this effect. Tricritical wings are predicted to appear 
in the phase diagram in a spatially varying magnetic field, but not in a homogeneous one. 
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Analogies between certain liquid-crystalphases (smec- 
tic, cholesteric, nematic, columnar, blue) [l| and ordered 
states of electrons in solids provide a suprising connec- 
tion between soft condensed matter and electronic sys- 
tems. For instance, stripe phases in high-temperature 
superconductors are analogous to smectics @, [1], helical 
magnets are analogous to cholesterics and electronic 
analogs of nematics can be realized via Pomeranchuk in- 
stabilities of the Fermi surface Q . Stripe and nematic or- 
der can occur in the spin channel as well as in the charge 
channel, and spin nematics provide magnetic analogs of 
non-s-wave superconductors that have been invoked to 
explain the enigmatic behavior of Sr3Ru207 @, 0] and 
the 'hidden order' in URU2S12 @|. 

For a theoretical description of phase transitions, Lan- 
dau theory is a standard tool [9j. Once an order pa- 
rameter (OP) <f> has been identified one constructs all 
terms in the free energy that are consistent with the gen- 
eral symmetry properties of the system. Replacing the 
OP field by a constant yields the most general mean- 
field theory for the phase transition. Fluctuations of the 
OP can also be considered, which leads to the Landau- 
Ginzburg- Wilson (LGW) framework that is amenable to 
an analysis by renormalization-group methods fioj ] . Lan- 
dau theory is applicable to both classical transitions that 
are driven by thermal fluctuations and tuned by chang- 
ing the temperature T, and quantum transitions at zero 
temperature that are driven by quantum fluctuations and 
tuned by a non-thermal parameter such as external pres- 
sure or composition [ill ]. 

In order to experimentally check the proposals con- 
cerning liquid-crystal analogs in electronic systems it is 
important to predict as many qualitative features of this 
exotic order as possible. To this end a comprehensive 
analysis of the Landau theory for an electronic nematic 
state has been developed, and the quantum phase tran- 
sition has been analyzed in analogy to Hertz's theory 
of s-wave quantum ferromagnets |12 |. see Ref. [HI and 
references therein. The theory has been developed for 
both the spin and the charge channel; here we will fo- 



cus on the spin channel and especially on the p-wave 
case. The OP for a non-s-wave ferromagnet is of the 
form {yJ} a {x)cj' l ab f{d x )yb} > {x)), with a 1 ' 2,3 the Pauli ma- 
trices, vb = ipi) is a fermionic spinor field, ip is its 
adjoint, x = (x,t) comprises the real space position x 
and the imaginary time variable r, / is a tensor- valued 
monomial function of the gradient operator, and (. . .) de- 
notes the quantum mechanical and thermodynamic av- 
erage. In the p-wave case the function / is linear in the 
gradient and the OP field thus carries a spin index i and 
an orbital index a (i = 1,2,3, a = 1,2,3). We will de- 
note it by Ni_ a (x) which in Landau becomes a number 
Ni ta . The most general Landau free energy up to quartic 
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tN?Ni + ui (KNtf + u 2 N?Ni NlN?, (1) 



with a summation convention implied. Depending on the 
relative values of the Landau parameters u\ and U2 there 
arc two distinct phases, the a- phase and the /3-phase fl3j . 
They are characterized by N^ a — n a Ni for the a-phase, 
with N a unit vector in spin space, and Ni 
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the /3-phase, with N a number and e^ 1 ' 2,3 * 1 three mutually 
orthogonal unit vectors. The Landau theory predicts a 
second-order transition into either phase at t = 0. The 
corresponding phase diagram is shown in Fig. [TJ For the 
classical transition at T > in 3-d fluctuations of the OP 
will modify the mean-field critical behavior predicted by 
Eq. (JlJ . For the corresponding quantum phase transition 
the dynamical critical exponent stabilizes the Gaussian 
fixed point in 3-d and the LGW theory that generalizes 
Eq. ([1]) predicts the transition to be second order with 
mean-field critical behavior (l3| . 

There are reasons to doubt the validity of the latter 
prediction. Despite its generality, Landau theory hinges 
on certain implicit assumptions. One is that an expan- 
sion of the free energy in powers of the OP is well be- 
haved, i.e., that the coefficients in a Taylor expansion are 
finite. The analogous statement within an LGW frame- 
work is that the field theory is local. This assumption can 
be violated if there arc soft modes other than the OP flue- 
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FIG. 1: Phase diagram for the Landau free energy given in Eq. 
|T]). The a and /3-phases are stable for U2 > 0, 3ui > — U2, 
and U2 < 0, Mi > — U2, respectively. The Landau theory 
predicts a second-order transition into either phase at t = 0. 



tuations that couple to the OP. In such a situation, the 
construction of a Landau theory entirely in terms of the 
OP amounts to integrating out these soft modes, which 
may lead to a non-local theory. Technically, a nonzero 
OP can give the additional soft mode(s) a mass, which 
implies that the free energy as a function of the OP can- 
not be analytic at <f> = 0. The resulting changes to the 
phase transition depend on the properties of the original 
Landau theory and on the sign of the leading nonana- 
lytic term. A classic example is the case of a type-I or 
weakly type-II BCS superconductor, where a simple Lan- 
dau theory predicts a second-order or continuous tran- 
sition. However, the coupling of the OP to the electro- 
magnetic vector potential leads to an effective free energy 
that contains a term of cubic order in the OP, which leads 
to a first-order phase transition [13]. The same is true 
for the nematic-to-smectic-A transition in liquid crystals, 
with the nematic Goldstone modes playing a role analo- 
gous to that of the vector potential. This phenomenon 
is known as a fluctuation-induced first-order transition in 
condensed matter physics, and as the Coleman- Weinberg 
mechanism for mass generation in particle physics [ToT ]. 
Another possibility is that the transition remains sec- 
ond order, but the coupling to the additional soft modes 
changes the universality class. This is believed to be the 
case deep in the type-II region (l6, 17|. 

For the quantum phase transition from a Fermi liq- 
uid to an electronic nematic state there are soft modes 
in addition to the OP fluctuations, namely, particle-hole 
excitations that are soft at T = and acquire a mass 
at nonzero T. The question is whether these soft modes 
couple to the OP in a way that invalidates simple Landau 
theory. In this Letter we show that they do for spin ne- 
matics (but not for charge nematics) and generically lead 
to a quantum phase transition that is fluctuation-induced 
first order. For simplicity, we will discuss the p-wave case 
in three dimensions (3-d) unless noted otherwise. 

We will first state and discuss our results, and then 
sketch their derivation. In a generalized mean-field the- 
ory analogous to Ref. [3 that neglects order-parameter 



fluctuations, but keeps the fermionic degrees of freedom 
and their coupling to the OP in a Gaussian approxima- 
tion one finds, for a 3-d system, an equation of state 
whose qualitative features are represented by 



>N 



tN + vN z \n(N 2 + T 2 )+uN 3 + o(A 3 ). (2) 



Here N is the number-valued OP for the a-phase (the re- 
sult for the /3-phase is structurally the same), u = U1+U2, 
and Hn is the field conjugate to N (more on this be- 
low), v > is a positive definite Landau coefficient that 
is given in terms of spin-triplet interaction amplitudes. 
o(x) denotes terms smaller than X ciS X ->■ 0. The Lan- 
dau parameters t and u are also (weakly) T-dependcnt, 
but we explicitly show only the T-dependence of the log- 
arithmic term since it cuts off a nonanalytic dependence 
on N . For 1 < d < 3 the corresponding nonanalyticity 
at T = is given by N d instead of A 3 In A. 

Equation @ predicts several features that are quali- 
tatively different from the Landau theory of Ref. |l3|. At 
T = 0, the term TV 3 In A is negative and larger than 
N 3 , which drives the transition first order below a tri- 
critical temperature T tc = exp(— u/2v). At T = there 
is a first-order transition at t = t± > that preempts 
the second-order transition predicted by Landau theory 
and leads to a discontinuous change of N from to 
Ni = exp[-(l+u/u)/2]. For h one finds t x = vNf/A. In 
the presence of a conjugate field, tricritical wings appear 
in the parameter space spanned by T, /ijv, and t, as is the 
case at any tricritical point 18]. As a result, we predict 
the phase diagram to have the same qualitative structure 
as the one observed for quantum ferromagnets 19|, 20| . 
see Fig. [2] There are surfaces ("wings") of first-order 
transitions that are bounded by second-order transitions 
and end in a pair of quantum critical points in the T = 
plane. The quantum critical exponents at these points 
are mean-field like. We choose as independent the static 
exponents /? and <5, and the dynamical exponent z, i.e., 
the scale dimension [T] c of T at criticality. We find 



0=1/2 



<5 = 3 



z = 3 



(3) 



These exponents are exact, since the conjugate field ren- 
ders massive the fermionic soft modes that, in the ab- 
sence of a field, lead to the first-order transition described 
above. /3 and S govern the dependence of N on t and hjy, 
respectively and z determines the temperature scaling of, 
e.g., the specific heat. The T-dependence of the OP at 
criticality, however, is not given by z for the same reason 
as in the s-wave case [U [22[ ■ It is determined by a fluc- 
tuation scale [T]fl uc = 9/(d+ 1) that yields the leading 
T-dependence for all d < 5. For d = 3 we find for the 
deviation SN of N from its value at the critical point 



6N(t c ,h c N ,T)ot-T 4 / 9 , 



(4) 



which also is exact. The continuous transitions at T > 0, 
on the other hand, will be modified by OP fluctuations. 
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FIG. 2: Schematic phase diagram in the space spanned by 
T, ftjv, and the non-thermal control parameter t. Shown are 
the region of p-wave order, the tricritical point (TCP), the 
two quantum critical points (QCP), and the various phase 
transition lines and surfaces. For /ijv = the transition is first 
order above the TCP and second order below. For h,N ^ 
tricritical wings connect the TCP and the two QCPs. 



In d = 2 there is no long-range order at T > 0. How- 
ever, at T — one still has a first-order quantum phase 
transition, with N x = (3v/4u) 2 , h =uNf/6. 

These predictions are all analogous to the properties of 
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quantum s-wave ferromagnets [n 
difference that allows to distinguish the two cases is the 
nature of the conjugate field. In the s-wave case, it is the 
physical magnetic field. In the p-wave case (as well as for 
all higher angular momenta) a homogeneous magnetic 
field does not couple to the OP. However, a conjugate 
field hjf is induced by a non-homogeneous magnetic field. 
For instance, a magnetic field h(x) with components hi 
in two independent directions induces a conjugate field 



h l ^(x) oc e ijk J dy hj(y) d a h k (x) 



(5) 



with V the system volume and a prefactor proportional to 
the p-wave interaction amplitude. A homogeneous field 
hpj results if hk{x) is a linear function of position. 

To summarize, our theory of low-T p-wave ferromag- 
nets or magnetic nematics predicts a first-order transition 
at asymptotically low T that is separated from a second- 
order transition at higher T by a tricritical point. In this 
respect it behaves the same way as an ordinary s-wave 
ferromagnet. However, in contrast to the s-wave case 
a p-wave ferromagnet is unaffected by a homogeneous 
magnetic field, whereas in a properly designed inhomo- 
geneous field it shows the same strikin g tr icritical-wing 
structure found in s-wave ferromagnets 13. 2pL 24j. 

A derivation of these results hinges on an effective the- 
ory of soft modes in clean electron systems that will be 



described in detail elsewhere [25[; here we just sketch 
some salient features. The building blocks are bilin- 
ear products ip ni (a?) ?/> n2 (y) that can be constrained to 
bosonic, quaternion- valued, fields Q nin2 (x,y). Here ri\ 
and ri2 label fermionic Matsubara frequencies u n — 
2irT(n + 1/2). This was pioneered by Wegner in the con- 
text of noninteracting electronic systems with quenched 
disorder [2f| ^gene ralized by Finkel'stein to interact- 
ing electrons I28L 1298 . and later elaborated on by the 
present authors [3(]j|. These theories were all formulated 
in terms of local density matrix fields Q ni n 2 i x i x )- Here 
we need to generalize to non-local phase-space variables 
Q ni n 2 y)- F° r a spherical Fermi surface we expand the 
Fourier transform of the Q in spherical harmonics Y™ , 



v p 



y/"(fi p )Q nin2 (p+fc/2,p-fc/2), 

(6) 

with ftp the solid angle for the wave vector p. For our 
purposes we need an effective theory for soft modes in 
clean systems. This is a harder problem than the dis- 
ordered one since in clean systems there are many more 
soft modes, namely, all moments of the phase-space exci- 
tations, not just the zeroth one. A Ward identity shows 
that all QiC^i 2 (k) with n\U2 < are soft. They are 
the Goldstone modes of a symmetry between retarded 
and advanced degrees of freedom that is spontaneously 
broken if the Fermi energy lies within a band; we de- 
note them by qll^^k) [31]. The role of the q in the 
present theory is analogous to that of the vector poten- 
tial in the case of the BCS transition in Ref. [lj. The 
various angular momentum channels, / = 0, 1, . . ., couple 
and all q-correlation functions scale as the inverse of the 
wave number which in turn scales as a frequency. In a 
schematic notation, leaving out all constant prefactors, 



J 7ll — 712 ,n>3 — 714 



\k\ + n r , 



(7) 



with f2„ = 2i:Tn a bosonic Matsubara frequency. This 
structure holds for both noninteracting and interacting 
electron systems, which reflects the fact that the sys- 
tem in the absence of magnetic order is a Fermi liquid. 
By integrating out the massive degrees of freedom in a 
conserving approximation one can construct an effective 
soft- mode theory in terms of the q only that is analogous 
to, albeit more complex than, the generalized nonlinear 
sigma-model for disordered interacting electrons [H, [3(| . 
The generalized mean-field theory discussed above can 
be derived analogously to the corresponding theory for 



the s-wave case 23 



32|. The leading coupling between 
the static, homogeneous OP N and the fermionic soft 
modes q takes the form N Tr [q(k) q^(k)] where Tr traces 
over the angular momentum, frequency, and spin degrees 
of freedom. To leading order in N and q one finds a 
fermionic part of the action with the schematic structure 



A q oc Tr 



\k\ + n + N) q(k) q^ {k)} 



(8) 
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Upon integrating out q one obtains a contribution to the 
equation of state that is schematically given by 



Nk 2 m 



, ^ (k 2 + n 2 n ) 2 (k + n n ) 2 + N 2 k 2 n 2 n 

(9) 

with a positive prefactor. I(N, T) appears in the equa- 
tion of state in addition to the standard Landau terms. 

Asymptotic analysis reveals that the singularity at 
N = 0, T = takes the form 



I(N,T = 0) oc N 



j const. - N^-V (l<d<3), 
1 const. + N 2 In N (d = 3), 



(10a) 



I(N, T)/N\ N _^ oc const. + T (1 < d < 3). (10b) 

For d — 3, and neglecting the analytic T-dependence 
of the term linear in N, the u-term Eq. © adequately 
represents the coupling of the soft fermionic modes to the 
OP (note the absence of a T 2 InT term in d = 3). 

We close with four remarks. First, the coupling of 
O(Nqq^) that generates the singular dependence on N 
is missing from the Hertz-type theory of Ref. [l3j. It 
amounts to taking into account fermionic loops in de- 
riving the LGW theory, whereas Hertz theory treats the 
fermions in a tree approximation. Second, our results do 
not carry over to charge nematics. In that case the OP 
does not give the q a mass and the Hertz-type LGW the- 
ory [l3T | is valid, consistent with a general argument first 
given in Ref. [33j. Third, the generalized mean-field the- 
ory still contains several approximations. The Gaussian 
approximation for the q is of no qualitative consequence; 
Fermi-liquid theory ensures that the exact propagators 
have the same scaling properties as Eq. 0. Higher-order 
coupling terms between N and the q are irrelevant in a 
renormalization-group sense and keeping only Eq. (j8j) suf- 
fices. This leaves the mean-field approximation for the 
OP. This becomes exact in a well-defined limit. N ^ 
gives the <7-propagator a mass, see Eq. ©, which defines 
a length scale A oc X/N, Together with the magnetic 
correlation length £ this defines a Ginzburg parameter 
k = A/£. For k — y (extreme type-I case) OP fluc- 
tuations are negligible and low-T transition inevitably 
is first order. For k > 0, and especially for k >> 1 
(type-II limit), the role of OP fluctuations must be exam- 
ined. This is a difficult problem that has been studied for 
classical transitions [3, H3] , and to some extent for the 
quantum s-wave ferromagnetic one [12] , but more work is 
needed on magnetic transitions in this limit. Fourth, all 
of our conclusions also hold for spin nematics in higher 
angular momentum channels, only the realization of the 
conjugate field changes. 
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